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Automatic Feynman diagram calculation with xloops 

-a short overview 
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Abstract 

xloops is a program package that calculates Feynman diagrams by using computer algebra 
systems. In this paper it is shown which problems to be solved by computer algebra arise 
during such calculations, and how this problems are handled in the framework of xloops. 

1 Introduction 



XLoops 1,1 
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Maple V loaded ! 

> read( '/home/Eruecher ALoop3/xloop3_l. l_aLpha/MapleVR3/loop3. ma 1 ) : 
Warning: new definition for 

:»»»»:>:>»»:»»:>:>»:»»»»»:*» This is XLOOPS Version 1.0 

by Lars Briicher. Johannes Franskowski, Dirk Kreimer 

©1993-1997 Johannes Gutenberg Universitat Mains, Germany 

bytes used=1006696 J 31100=851812, time=0.98 
bytes used=2020300 J alloc=1310480, time=2.42 
Model := 'SM' ; 

Model := SM 



Figure 1: The graphical user interface 



The increasing precision of experimental data in High energy Physics forces theorists to give 
very exact predictions for cross section of reactions between elementary particles. Calculating 
such cross sections efforts taking hundreds or thousands of Feynman diagrams into account. Every 
such diagram represents mathematically a 4-fold integral over internal momenta. Of course this 
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calculations cannot be done by hand any more. So the usage of a computer algebra system is 
vital. 



2 Structure of xloops 

Building a computer program doing automatic calculations of Feynman diagrams efforts knowledge 
in different sections of science. Beside physical understanding knowledge in mathematics and 
computer science is necessary to solve the arising problems. This together with just the number 
of problem to be solved makes it necessary to give xloops a structure with well defined interfaces 
between different programs. How this is realized in xloops can be seen in fig. |^. There can also be 
seen that xloops uses different programming languages for different tasks. First the graphical user 
interface (GUI) allows the user to input the process he is interested in, as easy and comfortable as 
possible. The main evaluation of the Feynman diagram is done by the MAPLE part of the program. 
Its major task is to produce either the full result for a Feynman diagram or a result that can be 
integrated numerically. Therefore xloops reduces the whole integral to sets of standard integrals, 
which can be solved in an algorithmic manner. To do so, xloops has to handle non-commutative 
objects like elements of the Dirac algebra and be able to evaluate traces of such objects. Finally 
these objects have to be exported to C++ code, so that they are ready for numerical analysis. 




generate diagrams ^^^,„ 

produces all diagrams to given external particles 

(MAPLE) 

evaluate diagrams E v.iGra P h 

parallel and orthogonal space splitting 



OneLoop 

oneloop 
functions 

(MAPLE) 



TwoLoop 

twoloop 
functions 

(MAPLE) 



numerical integration 



(C++/VEGAS) 



Figure 2: The program structure of xloops 



3 Physical and mathematical background 

To clarify the physical background of Feynman diagram evaluation the following section will give 
a brief overview. 

The quantized structure of micro cosmos is visualized by Feynman diagram containing closed 
loops of internal particles. But such a Feynman diagram is also a prescription for evaluating its 
contribution to the cross section or decay width of a physical process. Every line (propagator) 
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and every point (vertex) corresponds to a mathematical expression as described by the so called 
Feynman rules. Multiplying these terms gives the whole contribution of the diagram. Closed 
loops of particles require integration over the momenta of the particles forming the loop. So every 
Feynman diagram corresponds to a n-fold integral, with n depending on dimension in space time 
and number of closed loops. The principal analytical structure of these diagrams is: 



where I represents a Lorenz vector and represents its components. Expanding the numerator 
and writing all loop-momentum independent variables in front of each integral lead to an integral 
containing just products of components of loop momentum in the numerator and the so called 
Propagators Vi — [(I + qi) 2 + m|] in the denominator. These integrals are called tensor integrals 
in contradiction to scalar integrals, which have no loop momentum in the numerator. They can 
now be further reduced to a number of standard integrals with different techniques. 

The standard technique, called Passarino Veltman method, tries to reduce these integrals to 
scalar integrals. In an intermediate step a linear system of equations is built and has to be solved. 
This linear system of equations can be rather big and so the most time consuming part of tensor 
integral calculation. 

An alternative technique splits the Minkowski space time into the space spanned by the external 
momenta (parallel space) and its orthogonal complement (orthogonal space) . The integration over 
the orthogonal space can now be transformed to an rather simple integration over the surface of a 
hyper-sphere and a remaining radius integration. Beside integrations done with residue theorem 
the integral formula 



expresses the result in "/^.-functions. Between these functions exist a lot of relations, which allow 
to reduce all integrals to principally one basic 7?.-function. 

A third technique uses the fact that in parallel and orthogonal space every scalar product 
between an outer momentum qi and a loop momentum / just projects out one component of the 
loop momentum. Every such product occurs in one of the propagators. So every component of 
Ifj, from parallel space can be expressed in terms of progators and masses m,. The remaining 
components from orthogonal space can be expressed in terms of the last propagator from (^), 
which has no external momentum. So a recursive definition, which reduces every tensor integral 
to its corresponding scalar integral and number of integral with less propagators is gained. 

4 Exemplary CA problem: reduction of tensor integrals 

In this section some techniques used by xloops should be shown in more detail. As example serve to 
different methods for evaluating tensor integrals with xloops. The standard technique (Passarino 
Veltman method) is not used by xloops, as solving linear systems of equations with computer 
algebra methods is much slower than the algorithms discussed below. 

4.1 A fast technique: recurrence relations with ^-functions 

In the case of the one-loop two-point function the resulting ^-functions gained by integration 
fulfill the following recurrence relation: 



(61 +b2)Ht(h, &25 zi, Z2) = (h + b 2 + 1) TZ t (bi + 1, 6 2 ; zi, z 2 ) - tz\ 1Z t -i(bi, 6 2 ; zi, z 2 ) (3) 
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resulting in the basic 7£-function: 
U- e {-\ + £, 1; Zi, z 2 ) 
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This can of course be directly translated into the following MAPLE-code: 
R2 := proc(ind,bl,b2,zl,z2) local erg; 



-l+2e' 



# *** Parameter increasing *** 
elif bl<0 then 
erg := ( (2*bl+2*b2+2*ind) * R2(ind,bl+1 ,b2,zl ,z2) 

- 2*ind*zl* R2(ind-l,bl+l,b2,zl,z2) )/(2*bl+2*b2) ; 



elif ind=0 then 
erg := z2~eps/ (2*zl~ (2*eps) ) * 

( (1+Sqrt (l-zl/z2) ) * (1-Sqrt (l-zl/z2) ) * (2*eps) 
+(l-Sqrt(l-zl/z2))*(l+Sqrt(l-zl/z2))~(2*eps) ) ; 
RETURN ( erg ); 

fi; 

end: 

The recurrence relations are for most cases twice as fast as solving the linear system of equa- 
tions, where as intermediate step matrices have to inverted. 



4.2 Even faster technique: cancellation of Propagators 

Parallel and orthogonal space splitting lead for the one-loop two-point function to the following 
relation (here just for parallel space component of loop momentum): 

h = ^-[v 1 -r 2 + c] (4) 

in terms of integrals this means: 

f JL = A_ 

J ViV 2 2<z|| 

with a integration independent constant C.So this equation reduces the exponent of loop momen- 
tum of a tensor integral by one and produces additional, simpler integrals. Applying this formula 
n times ends with the scalar two-point function. This can also be written in a iterative algorithm: 

FastTensor2Pt := proc (pO ,pl ,q0 ,m_l ,m_2 , eps) 

local Cl,C2,DM,i0,il,i0s,cf 1 , cf 2 , cf 3,cf 4, cf 5, cf 6,ml ,m2,temp0 ; 

global rho; 
if not type (pi .even) then RETURN ( ); fi; 
cf4 := 0; 

for iO from 1 to pO do 
iOs := trunc((i0-l)/2) ; 
cf3 := 0; 

for il from to iOs do 
cf 1 : = binomial (iO-1 , 2*il) * (2*q0) " (2*il-p0) *Two20nePtFactor (2*il , pi , eps) ; 
cf2 := (Cl)-(iO-l-2*il)*Tadpole(m2,l+il+pl/2-eps) 

+(-1) ~iO*C2~ (i0-l-2*il) *Tadpole (ml , l+il+pl/2-eps) ; 



f 1 f 1 +c f JL 

J V 2 J Vi J Pi 
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cf3 := cf3 + cfl*cf2; 
od; 

cf4 := cf4 + binomial (pO, iO) * (-C1) ~ (pO-iO) *cf 3; 
od; 

cf6 := cf4/(l-eps) ; 

tempO := (-Cl/2/qO) ~pO*Small_2Pt (pi ,q0~2 ,ml ,m2, eps)+cf 6 ; 

tempO := evaKsubs (Cl=qO~2-m_l+m_2,C2=qO~2-m_2+m_l ,ml=m_l ,m2=m_2,temp0) ) ; 
RETURN ( eval(tempO) ); 
end: 

Depending on the dimension of parallel space this technique is even a bit faster than the above 
described. 

5 Conclusion and outlook 

In high energy physics different techniques of symbolic calculation are used. Some of them used 
by xloops have been pointed out here in detail. Moreover xloops uses the following elements and 
algorithm of computer algebra systems: 

• user defined recursive functions 

• power series expansion 

• fast list evaluation 

• algebraic term collection 

The hole package is still a developing project. It can be obtained from its homepage on the WWW: 



http : //wwwthep . physik . uni-mainz . de/^xloops 
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